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Abstract 


This  proposal  addresses  two  related  aspects  of  the  acoustic  reverberation  problem, 
namely  (1)  seismic-acoustic  properties  of  ocean  bottom/subbottom  sediments,  and  (2)  the 
short-range  scattering  of  acoustic/elastic  waves  at  the  rough  and  heterogeneous 
bottom/subbottom  interface  zone.  The  primary  objective  of  the  proposed  research  is  to  quantify 
the  importance  of  bottom/subbottom  scattering  of  compressional  and  converted  shear  waves. 
The  approach  is  to  numerically  simulate  short-range  scattering  using  realistic  models  for  the 
geometry,  scales,  and  vertical  and  lateral  variations  in  material  properties  of  the  ocean  bottom 
and  especially  subbottom.  To  treat  this  problem  realistically,  a  numerical  model  for  constrained 
sedimentological  deposition  of  sands  and  clays  over  basement  rock  is  used  to  describe  the  three- 
dimensional  geometry  of  sediment  beds  and  laminae  over  basement  rock.  Models  for  rock 
properties  are  then  used  to  determine  the  seismic  velocity,  attenuation,  density,  porosity,  and 
permeability  for  the  suspended  and  loosely  consolidated  granular  sediment  materials  and 
basement  rocks.  This  geological  model  is  then  used  in  three-dimensional  seismic-acoustic 
simulators  to  predict  the  acoustic  fields  scattered  by  the  rocks  at  and  below  the  interface.  Two 
preliminary  finding  of  results  by  us  to  date  are  noteworthy:  First,  acoustic  scattering  is 
completely  inadequate  for  describing  the  interaction,  particularly  near  grazing  angles  where  the 
elastic  boundary  conditions  play  a  critical  role  in  establishing  the  energy  balance.  Second, 
bathymetrically  "smooth"  surfaces  may  be  acoustically  very  "rough",  providing  a  mechanism 
for  the  "anomalous"  scattering  regions  observed  in  reverberation  data  and  the  importance  of 
subsurface  scattering  in  the  reverberation  problem.  Our  simulation  codes  are  implemented  on  a 
dedicated  in-house  supercomputer,  thus  providing  a  powerful  and  cost  effective  computer 
platform  for  performing  the  proposed  research. 
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1.0  Science  Statement  and  Research  Goals 


For  purposes  of  describing  the  reverberant  acoustic  wavefield  in  the  ocean  column,  the 
seismic-acoustic  scattering  phenomena  at  the  bottom/subbottom  interface  can  be  conveniently 
divided  into  two  range-scale  dependent  regimes:  (1)  short-range  (several  tens  of  wavelengths) 
interactions  at  the  bottom/subbottom  interface,  e.g.,  scattering  zones,  tmd  (2)  long-range 
hundreds  of  kilometers)  interactions  among  spatially  distributed  scattering  zones.  In  this 
proposal  we  address  the  short-range  scattering  phenomena  using  2D  and  3D  numerical 
simulators  of  scattering  from  realistic  geological  models  for  ocean  bottom/subbottom  sediments. 

Preliminary  research  in  this  area  paformed  in  our  group  (Samec  and  Marion,  1989)  has 
revealed  two  important  physical  phenomena  related  to  sea  floor  scattering.  First,  we  find  that  a 
heterogeneous  shear  transition  zone,  where  the  shear  velocity  increases  slowly  from  zero  to 
some  finite  value,  must  exist  between  the  water  and  the  compacted  sediments.  Acoustic  waves 
penetrating  the  bottom  are  continuously  scattered  in  diis  zone  with  some  enCTgy  being  converted 
to  downgoing  forward  scattered  shear  waves.  In  a  similar  way,  distributed  scattering  of 
upgoing  compressional  and  shear  waves,  e.g.,  reflected  from  basement  rocks,  generate  upgoing 
compressional  energy  which  reenters  the  ocean  column  and  contributes  to  the  reverberant 
acoustic  field.  It  is  therefore  critically  important  to  consider  the  full  elastic  scattering  which 
includes  shear  waves  if  the  scattered  energy  balance  is  to  be  quantitatively  determined. 

The  second  important  observation  from  our  preliminary  study  is  that  a  bathymetrically 
"smooth"  surface,  as  generated  by  the  SedSim  depositional  model  described  below,  can  appear 
very  "rough"  to  the  illuminating  acoustic  wavefield.  This  effect  depends  on  the  history  of 
deposition  and  compaction  of  the  sediments,  the  frequency  of  the  acoustic  waves,  the  angle  of 
incidence  of  the  acoustic  wavefield,  and  of  course,  the  material  properties  and  the  spatial 
variability  of  the  sediments.  The  question  arises  as  to  whether  this  effect  is  partially  responsible 
or  important  for  explaining  the  anomalous  zones  called  "phantoms"  which  are  often  observed  in 
real  reverberation  data. 


The  observations  summarized  above  were  made  using  2D  geological  elastic  models.  The 
questions  raised  in  the  SRP  announcement  and  subsequent  planning  meetings  can  only  be 
answered  with  realistic  geological  and  seismic-acoustic  models  -  three-dimensional  anelastic 
scattering  from  realistic  geological  models  not  only  with  rough  water  bottom  but  also  with 
subbottom  heterogeneities. 


The  primary  goal  of  the  proposed  research  is  to  quantitatively  describe  the  importance  of 
short-range  bottom/subbottom  scattering  of  compressional  and  shear  waves^n  the 
compressional  wave  energy  reentering  the  ocean  column.  This  research  has distinct 
components:  ^ 


(1)  Realistic  geological  modeling  to  describe  the  rough  interface  and  subbottom 
volume  heterogeneities  and  rock  physics  modeling  to  describe  material  properties  of  the 
sediments. 


(2)  Seismic-acoustic  modeling  to  identify  the  role  of  subbottom  scattering  on  the 
acoustic  energy  balance,  closely  examining  the  influence  of  shear  waves  and  attenuation. 

(3)  Statistical  analysis  of  the  short-range  scattered  fields  for  use  in  long-range 
reverberation  models. 
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2.0  Research  Approach  and  Preliminary  Results 

2.1  Technical  Approach 

A  sienificant  effort  is  already  underway  in  our  research  groups  in  areas  of  numerical 
modeling  of  seismic  scattering  in  heterogeneous  media.  A  variety  of  tectaiques  with  worlang 
computer  codes  are  in  place.  These  include  3D  ray  ttacers  with  pmbohc  approximations 
gaussian  beams,  and  full  wave  solvers  using  finite  ^fferences.fmite  dements, 
techniques.  Research  is  being  started  on  boundary  eleinent  methods.  This  ongoing  research 
has  mostly  been  directed  at  and  funded  by  exploration  and  development  applications  for  oil  ^ 
gas.  We  propose  to  adapt  the  existing  2D  and  3D  full  wave  seismic  modeling  code^e.g.,  re, 
FD  and  Hybnds,  to  the  problem  of  short-range  scattering  firom  the  ocean  bottom  mterface.  The 
oon^onvrocts,  is  akeady  underway  under  FY89  ONR  seed  funding.  By  short-r^  we 
mean  3D  dimensions  of  the  order  of  50x50x50  wavelengths.  A  summary  of  the  computer 
algorithms  available  for  this  project  is  given  in  Table  1. 

Similarly,  we  have  existing  models  for  sediment  deposition  and  rock  prop^es  which 
together  generate  realistic  geological  models  for  the  bottom/subbottom  interface.  By  realisnc, 
we  mean  that  the  scales  of  heterogeneities  and  material  properties  are  consistent  with  geological 
constraints  and  experimental  observations. 

Both  geological  and  seismic-acoustic  models  are  run  on  our  in-house  Ardent  Titan-O 
supercomputer.  This  machine  is  wholly  owned  by  the  Department  of  Geophysics  thus 
providing  a  powerful  in-house  platform  for  the  numencal  work  with  rapid  access  to  resul  s. 
Our  version  of  the  Ardent  has  two  CPU’s  each  with  vector  P^ocf  sors  which  can  sustsun 
Mflops.  Estimates  for  the  CPU  time  needed  for  problems  *at  will  1^  address^  in  this  project 
are  given  in  Table  1.  The  preliminary  results  summanzed  below  m  Section  2.1  md  in  1  able  l 
indicate  that  our  approach  is  an  accurate  way  of  examining  both  2D  and  3D  scattenng  effects. 

There  are  two  major  issues  addressed  in  the  numerical  approach  used:  First,  an  accurate 
representation  of  physical  properties  of  the  geological  model  -  such  as  anisotropy  and 
attenuation  of  unconsolidated  sediments,  and  second,  an  adequate  way  of  modeling  the  seismic 
response  to  geometry  of  the  water  sediment  interface  and  the  small  scale  volumetric 

heterogeneities. 

Arniratelv  MoHp.ling  the  Sediments.  To  describe  the  scattering  interaction,  one  must 
parameterize  the  properties  of  the  ocean  bottom  and  subbottom  with  elastic  m^uli,  attenuation 
factors  (6),  density,  and  a  geometric  description  of  the  heterogeneities  including  the  rough 
interface  SedSim  (Harbough,  Dept  of  Applied  Earth  Sciences)  is  a  numencd  m^el  for 
constrained  sedimentological  deposition  of  sand-clay  concentrations  representative  of  yoimg 
marine  sediments.  Multiple  depositional  episodes  can  be  simulated  to  form  the  3D  geometty, 
i  e ,  structure,  of  beds,  laminae,  and  ponds.  The  scales  of  heterogeneities  generated  are  sel^ 
similar  so  that  the  models  are  dimensionless.  Also,  for  several  years  ^d  with  continuous  ONR 
support,  our  Rock  Physics  Group  has  been  developing  seisimc  velocity  ^d  Q  models  for 
sedimentary  mixtures  of  sand  and  clay.  Based  on  both  field  and  laboratory  results,  these 
models  describe  relationships  between  seismic  velocity,  density,  attenuation,  pomsity, 
permeability,  mineralogy,  and  confining  pressure  for  suspended  and  loosely  consohda  e 
granular  materials.  Together,  the  depositional  model  and  the  seismic  sedinient  propemes 
models  provide  for  the  first  time  the  necessary  tools  for  creating  geologically  realistic 
sedimentary  models  for  the  ocean  bottom/subbottom  interface. 

Arniiratelv  Mndp.ling  the  Acoustic  Scattering.  Both  finite  difference  and  Fourier 
algorithms  provide  us  with  fast  and  computer  efficient  methods  to  evaluate  spatial  operators  tor 
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material  properties  that  are  adequately  sampled  on  regular  grids  (Kosloff  1973).  The  Fourier 
algorithm  tends  to  provide  a  more  accurate  evaluation  of  spanal  operators  for  models  wth 
rapidly  varying  properties  as  in  the  case  of  the  SedSim  model  investigated  to  date. 
Nevei4elesi  both  algorithms  deliver  excellent  accmcy  for  four  or  more  gnd  ^ints  per 
wavelength.  In  our  existing  codes,  both  finite  (Hfference  and  Fourier  algonthms  are 
implemented  with  a  central  difference  method  for  time  integration. 

Modeling  of  attenuation  is  performed  using  a  pseudo  spectral  algorithm  (Kosloff  1988, 
Samec  1987,  Carcione  1987).  The  attenuation  mechanisms  are  modeled  by  a  set  of  relaxation 
mechanisms  (two  or  more)  therefore  providing  a  means  to  shape  the  variations  of  the  quality 
factor  (Q)  with  jfrequency  according  to  experimental  results,  and  still  yieldmg  a  causal 

attenuation  mechanism. 

Numerical  modeling  of  regularly  sampled  data  with  the  algorithms  inentioned  above 
results  in  well  known  artifacts  for  curved  interfaces.  Such  artifacts  are  numenc^ly  overwme 
by  using  finite  element  techniques  which  allow  an  accurate  description  of  both  interface 
geometnes,  as  well  as  an  explicit  formulation  of  boundary  conditions.  To  properly  h^dle  both 
the  geometrical  features  of  the  sea  floor  and  the  scales  of  heterogeneities,  an  exphcu-imphci 
finite  element  (Hughes  1983)  algorithm  has  been  implemented  that  allow^he  h^dling  of  a 
large  number  of  elements  (35,000)  in  acceptable  y^dent  cornputetion  time.  perform^ce  of 
this  algorithm  is  further  improved  by  the  use  of  an  hybnd  Finite-elemen^nite-differenw 
algorithm  Aat  both  allows  a  careful  description  of  the  geometncal  features  of  the  oceM  bottom 
as  well  as  a  fast  and  accurate  numerical  modeling  of  the  subbottom  volume  heterogeneities. 

2.2  Preliminary  Results 

(The  results  shown  here  were  obtained  in  part  in  early  July  under  the  auspices  ^ 
anticipated  FY89  0NR  funding.)  The  SedSim  depositional  simulator  was  run  to  create  a  3D 
geological  model  for  the  sediments.  In  Fig.  1  left,  we  illustrate  three  2D  slices  firorn  the  3D 
SedSim  model.  The  water  column  is  shown  in  white  overlying  the  sedirnents  built  up  on 
basement  rock.  The  rock  properties  model  was  then  applied  to  cornp^t  the  segments  and 
assign  values  for  density,' compressional  and  shear  wave  velocities.  In  Fig.  1  nght,  we  show 
the  elastic  properties  for  SedSim  slice  A.  Note  the  fuzzy  shear  trpsition  zone  ^  nght)  near  Ae 
water-sediments  interface.  The  compressional  wave  reflectivity  is  shewn  in  Fig.  1C  ng  . 
Reflectivity  spectra  through  this  model  are  shown  in  Fig.  2  to  illustrate  the  self-similar  power 
law  distribution  of  reflectors,  seen  by  acoustic  waves,  in  the  heterogeneities  generated  by 

SedSim. 

In  the  next  two  figures  we  display  the  scattered  wavefield  for  simulated  experiments  run 
over  the  SedSim  model.  Fig.  3  is  acoustic;  Fig.  4  is  elastic.  The  experiment  has  an  array  of 
detectors  collinear  with  an  isotropic  source  in  the  water  above  the  sediments.  In  these  and 
subsequent  seismograms,  we  plot  the  detected  displacement  in  water  rather  than  the  pre^e; 
this  results  from  our  seismic  applications,  a  result  which  will  be  changed  as  part  of  our  r  loy 
research.  In  both  figures  3  and  4,  the  direct  wave  is  present  on  the  honzontal  detectOTS  only  -  as 
it  should  be  because  the  displacement  detectors  are  horizontally  collinear  with  with  the  source. 
The  elastic  scattering  wavefield  (Fig.  4)  is  considerably  more  complicated  than  the  acoustic 
wavefield  (Fig.  3),  particularly  in  the  later  arriving  events  where  the  sh®^  waves  contnbute 
most.  By  subtracting  the  acoustic  and  elastic  seismograms,  we  can  highlight  th^ner^which 
contributes  to  the  reverberation  which  is  neglected  by  the  acoustic  formulation,  ^s  difterence 
seismogram  is  shown  in  Fig.  5.  A  large  part  of  the  direct  wave  is  annihilated  by  the 
subtraction,  but  near  grazing  angles  the  acoustic  and  elastic  fields  differ  most  significantly. 

Clearly  the  elastic  formulation  is  needed  to  accurately  model  the  near  gazing  angle 
phenomena  so  important  to  the  reverberation  problem.  This  simple  subtraction  illustrates  t  e 
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importance  of  treating  elastic  propaties  and  qualitatively  de^bes  the  particular  role  play^  by 
shear  waves  generated  within  the  sediments.  Fu^ermore,  tecause  the  convert^  shear 
energy  reenteiing  the  ocean  tends  to  be  less  coherent  than  other  amvals,  it  is  likely  to  conmbute 
significantly  to  the  reverberations.  These  m^els  were  run  without  attenuation,  which  will 
furtfiCT  affect  Ae  relative  magnitude  of  shear  contribution. 


In  Fig.  6,  we  display  a  constant-offset  seismogram  -  where  the  source  a^detector  are 
set  a  fixed  horizontal  distance  apart  and  moved  laterally  across  the  sediments,  ms  ge^etry 
highlights  the  effects  of  lateral  heterogeneities  as  sensed  acoustically.  In  Fig.  6A,  the  offset  is 
oidy  50  meters  thus  giving  a  near  vertical  incidence  backscattered  unag®  of  the  sediinent 
structure.  Note  the  incoherence  of  the  returns.  In  Fig.  6B,  the  honzon^  offtet  is  350  meters, 
thus  approaching  grazing  incidence  forward  scattering.  In  this  case  the  reflections  from  the 
sedi^nts  appear  more  coherent,  an  effect  due  to  lateral  homogenization  introduced  by  the  near 
grazing  propagation. 


Finally  in  Figure  7,  we  show  a  result  for  attenuation  and  anisotropy  although  modeled 
here  only  in  a  homogeneous  medium.  This  result  is  included  simply  to  illustrate  ^e  kind  of 
effects  which  give  rise  to  Q  anisotropy  and  to  document  the  ability  of  our  codes  to  mclude  Q  and 
anisotropy  both  of  which  will  be  used  in  the  ocean  bottom  research. 
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Figure  1:  Sedsim  model  :  Displayed  on  the  left  are  three  slices  extracted 
from  the  Sedsim  three  dimensionnal  sediment  modd.  From  A  to  C,  the 
slices  are  extracted  further  from  the  sediment  input  location  (further 
from  the  coast  line),  which  explains  the  thinning  of  the  sediment  layer. 
Displayed  on  the  right  are  the  compressionaJ  velocity  (A),  the  shear  ve¬ 
locity  (B),  and  the  reflection  coefficient  (C)  inferred  from  sediment  bulk 
composition. 
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Reflectivity  Spectrum 


Figure  2:  Frequency  distribution  of  the  the  reflectors  as  seen  in  Fig.  Ic:  This 
graph  presents  the  power  spectrum  of  the  reflectivity  as  a  function  of  the 
two  way  travel  time.  The  straight  line  underlines  the  proportionality  of 
the  power  spectrum  to  ®  which  is  consistent  with  the  data  inferred 
from  seismic  and  well  log  data  in  Gulf  coast  sediments. 
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Figure  3:  Common  shot  gather  for  the  acoustic  modeling  :  a)  Horizontal 
component  of  the  displacements,  b)  Vertical  component  of  the  displace¬ 
ments. 
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Figure  4:  Common  shot  gather  for  the  elastic  modeling  :  a)  Horizontal  com¬ 
ponent  of  the  displacements,  b)  Vertical  component  of  the  displacements. 
Note  here  the  importance  of  the  shear  waves  related  events. 
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Figure  5:  Straight  amplitude  difference  between  the  models  presented  in  Fig. 
3,  and  Fig.  4. 
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Figure  6:  Constant  Offset  sections  :  a)  50  meters  offset,  b)  350  meters  offset. 
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3.0  Proposed  Activities  and  Schedule  For  FY90  and  FY91 

We  include  the  activities  proposed  for  the  FY89  research  program  in  order  to  show  the 
continuity  with  proposed  FY90  and  FY91  research.  A  time  line  schedule  of  activities  is  shown 
in  Figure  8. 

FY89:  Activities  during  FY89  is  largely  investigatory,  aimed  at  identifying  the  important 
short-range  scattering  questions  for  long-term  study. 

1 .  Use  2D  SedSim  to  generate  examples  of  the  lateral  and  depth  variations  representative  of 
the  scales  of  heterogeneities  found  in  young  marine  sediments. 

2.  Use  existing  rock  physics  models  to  include  compaction  in  the  SedSim  sediments  and 
assign  seismic  values  for  compressional  and  shear  velocities,  density,  porosity,  and 
permeability. 

3.  Assess  the  numerical  accuracy  of  the  2D  seismic  code  for  modeling  the  shear  transition 
zone  (shear  velocities  approaching  zero)  at  the  ocean  bottom/subbottom  interface. 

4.  Initiate  study  of  the  role  of  converted  shear  waves  and  attenuation  on  the  reradiation  of 
subbottom  scattered  energy  back  into  the  ocean  column. 

5.  Determine,  using  SedSim  models,  the  parameters  needed  to  effectively  chmeterize  the 
scattering,  e.g.,  velocity,  attenuation,  density,  or  impedance;  determine  the  media  parameters 
which  control  the  magnitude  of  the  water  /sediment  wave  interaction. 

6.  Examine  the  scales  of  heterogeneity  in  2D  and  3D  which  dominant  the  short-range 
scattering  interaction. 


FY90:  We  suggest  that  the  2D  seismic-acoustic  models  will  not  provide  good  quantitative 
representations  of  the  scattering  primarily  because  3D  geological  models  are  often 
microscopically  anisotropic  often  with  strong  lineations  perpendicular  to  the  direction  of 
depositional  flow.  Such  anisotropic  features  are  well  represented  in  the  SedSim  models.  This 
year’s  activities  include: 

1.  Expand  2D  elastic  simulations  started  in  FY89  to  3D  anelastic  simxilations.  This  requires 
improved  development  of  attenuation  models  for  sediments  using  lab  and  field  data  for 
control.  Refinement  and  optimization  of  the  seismic-acoustic  simulators  will  continue  with 
specific  focus  on  the  hybrid  finite  element-  finite  difference  algorithm. 

2.  Begin  modeling  the  canonical  revCTberation  experiment  with  actual  source  and  detector 
geometries,  frequencies,  and  beamwidths  for  source  and  detector  arrays.  Analysis  tools  for 
beamforming  point  source  and  point  detector  simulation  data  will  be  be  developed. 

3.  Using  the  footprint  formed  by  the  source  and  detector  array  as  a  guide  for  scale,  we  will 
begin  the  development  of  techniques  for  partitioning  the  scatter^  wavefield  into  deterministic 
and  stochastic  components.  An  expect^  output  of  this  activity  is  a  statistical  description  of 
the  differential  scattering  cross-sections  for  elastic  scattering  off  the  sediments.  Results  will 
be  compared  with  analytic  treatments  of  scattering  by  rough  surfaces  and  random  media. 
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FY91:  To  this  point,  we  have  developed  the  basic  simulation  ^d  analysis  tools  for  short- 
range  scattering.  Research  will  continue  on  improving  and  refining  these  tools.  A  specific 
activity  for  FY91  will  be  to  simulate  the  natural  lab. 

1 .  Complete  development  of  sediment  properties  for  velocity,  attenuation,  density,  porosity, 
and  permeability  and  the  respective  anisotropies  in  sedimentary  mixtures. 

2.  Finalize  the  seismic-acoustic  modeling  suite  for  3D  anelastic  media. 

3.  Develop  methods  for  embedding  the  short-range  scattering  results  into  the  long-range 
reverberation  models.  This  will  be  done  with  outside  researchers  yet  to  be  identified  who  are 
working  on  the  long-range  reverb  problem. 

4.  Create  and  run  specific  3D  deterministic  simulations  for  a  site  within  the  natural 
laboratory. 


4.0  Deliverables 

Several  ONR  funded  research  groups  will  be  addressing  the  short-range  scattering 
problem.  It  is  impcKtant  that  each  group  compare  results  so  that  a  final  assessment  can  be  made 
for  the  strengths  and  weaknesses  of  the  many  different  modeling  schemes  iinplemented.  To  that 
end,  we  propose  to  make  available  to  the  scattering  community  a  3D  SedSim  generated  model 
with  complete  anelastic  properties  for  the  sediments,  including  realistic  compressional  and  shear 
wave  velocities,  density,  permeability,  porosity,  and  eventutdly  compressional  and  shear  wave 
attenuation  factors.  This  model  with  have  both  topographic  roughness  at  the  sediments’ 
interface  and  volumetric  distributed  heterogeneities. 

The  major  deliverable  shall  be  a  suite  of  computer  codes  capable  of  modeling  the  short- 
range  scattering  of  seismic-acoustic  waves  for  three-dimensional  anelastic  media.  In  addition  to 
these  basic  wave  equation  solvers,  we  will  have  methods  and  codes  for  modeling  the  actual 
source  and  detector  georrietries  used  in  the  reverb  experiment,  that  is,  beamfotmed  source  and 
detectors  arrays.  Finally,  we  will  provide  methods  and  codes  for  analyzing  the  simulated 
results;  these  will  include  codes  for  calculating  differential  scattering  cross-sections  and 
statistical  properties  of  simulated  scattering  seismograms  and  acoustograms. 
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Schedule  of  Activities 
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Figure  8 


BUDGET 


10/1/89  -  9/30/90  10/1/90  -  9/30/91 


A.  SALARIES  AND  WAGES 


Amos  Nur,  P.I.,  .5  months, 
calendar  year 

$3,900 

$4,000 

Jerry  Harris,  Co-P.I.,  .5  months. 

calendar  year 

Postdoc  research  fellow,  5  mo 

's. 

3,900 

4,000 

100%  time 

Graduate  students  (one  first 

year. 

13,000 

13,000 

two  second  year)  50%  academic 
100%  summer 

yr.. 

14,000 

28.000 

TOTAL  SALARIES  AND  WAGES 

$34,800 

$49,000 

B. 

BENEFITS 

9/1/89  -  8/31/90  @  27.8% 
9/1/90  -  8/31/91  @  28.6% 
9/1/91  -  8/31/92  @  29.3% 

$9,699 

$14,043 

TOTAL  SALARIES  AND  BENEFITS 

$44,499 

$63,043 

C. 

TRAVEL 

Two  trips:  east  coast  &  west 

coast 

$2,000 

$3,000 

D. 

EXPENDABLES 

Computing  costs 

200  hours  @  $50/hour 

300  hours  @  $50/hour 

$10,000 

$15,000 

Publication  costs 

25  pages  @  $100/page 

50  pages  @  $100/page 

2,500 

5,000 

TOTAL  EXPENDABLES 

$12,500 

$20,000 

E. 

TOTAL  DIRECT  COSTS 

$58,999 

$86,043 

F, 

TOTAL  INDIRECT  COSTS 

7A%  of  MTDC 

$43,659 

$63,672 

G. 

TOTAL  COSTS 

$102,658 

$149,715 

H. 

TOTAL  TWO-YEAR  COSTS 

$252,373 
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STANFORD  UNIVERSITY 

STANFORD,  CALIFORNIA  94305-2215 


DEPARTMENT  OF  GEOPHYSICS 
School  of  Earth  Sciences 


February  20, 1990 


Mr.  Marvin  Blizard 
bct#l,  Rm633 
Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217-5000 


Ref:  ONR  grant  with  a  start  date  of  01-01-90  for  two  years  @  $252,373. 

Proposal  #  8942—0169  Procurement  #  4258043 — 01. 

Title:  Ocean  bottom/subbottom  seismic-acoustic  scattering  with  realistic  rediment 
properties. 

Principal  Investigators:  Jerry  M.  Harris/Amos  Nur 
Department:  Geophysics 
Dear  Mr.  Blizard: 

As  per  our  conversation  today,  this  is  a  request  for  a  budget  modification.  Attached 
please  find  a  modified  budget  for  the  above  mentioned  grant.  The  modification 
reflects  a  more  realistic  use  of  our  proposed  dollar  expenditure.  This  request  to 
rebudget  the  salary  portion  of  the  grant  will  not  effect  the  total  cost  of  the  project. 

Sincerely, 

u 

Margaret  Muir 

AA  to  Amos  Nur,  Chairman  and  Professor 
Geophysics  Department. 

end:  revised  budget. 

cc:  Stephanie  Mitchell,  Contracts  officer,  ONR. 

Lilly  Ryan-Culclager,  Jerry  M.  Harris 


MODinCATION  February  1, 1990 


A. 


B. 


C. 

D. 


E. 


F. 


ONR  BUDGET: 

01/01/90. 

01/01/91. 

S ALARffiS  AND  WAGES 

Jerry  Harris,  PL,  .05  months/calender 

Amos  Nur,  CO  PL,  .05  months/calender 

3,900.00 

3,900.00 

4,000.00 

4,000.00 

Postdoc  research  fellow,  5  mths  @  100% 

Graduate  students,  50%  academic  yr. 

(2  first  year,  3  second  year)  100%  summer. 

13,000.00 

22,000.00 

14,000.00 

39,000.00 

Secretary 

750.00 

800.00 

43,550.00 

61,800.00 

BENERTS  @ 

1st  year:  8  mths  @  27.6%  4  mths  @  28.3% 

2  year:  8  mths  @  28.3%  4  mths  @  29.6% 

12,119.97 

17,242.20 

TOTAL  SALARIES  AND  BENEFITS 

55,669.97 

79,042.20 

TRAVEL  (National  meetings). 

2,000.00 

3,000.00 

SERVICES  &  EXPENDABLES. 

Office/lab  supplies,  xeroxing,  mailing. 

881.33 

2,067.35 

TOTAL  DIRECT  COSTS 

55,669.97 

84,109.55 

INDIRECT  COSTS® 

1st.  year:  8  mths  @  74%,  4  mths  @  78% 

44,106.69 

65,605.45 

2nd.  year:  @  78%,. 

YEARLY  TOTALS: 

$102,657.99 

$149,715.00 

TOTAL  PROJECT  COST:  $252,373. 


6.0  Other  Anticipated  Support 

There  are  several  active  on-going  projects  within  the  Department  of  Geophysics  in  rock 
physics  and  seismic  modeling.  In  general,  these  projects  encompass  the  development  of 
seismic  modeling  codes  in  combination  with  research  in  seismic  inversion,  data  processing,  and 
data  acquisition.  Most  are  funded  by  the  oil  and  gas  industry.  In  addition,  Professor  Nur  has 
existing  DOE  and  ONR  funding  for  rock  properties  research.  Professor  Harris  has  support 
($940k  with  Nur)  from  the  Gas  Research  Institute  for  cross-well  seismic  and  anticipates  support 
($20k)  from  British  Petroleum  for  interactive  seismic  imaging  research. 


_ 

Project/grant 

Funding  Agency 

Period 

Funding 

PI  time 

Academic  yr. 

Physical  Properties 

CNR 

10/1/89-9/30/90 

75k 

8% 

of  Marine  Sed. 

Rocks  with  fiuids 

DC3E 

10/1/89-9/31/96 

175k 

8% 

Block  tectonics 

NASA 

3/1/89-2/28/90 

105k 

8% 

Rock  Physics  & 

GRi 

10/1/89-10/31/90 

115k 

8% 

Cross-weii  imaging 

(with  Harris) 

1/1/89-12/31/89 

1  946k  1 

1  8% 

7.0  Qualifications 

7.1  Facilities  and  Equipment 

Computing  facilities:  We  have  in  operation  an  ARDENT  Titan  graphics  supercomputer, 
which  provides  the  intense  numerical  computing  needs  for  the  project.  We  also  have  three  4 
SUN  workstations  and  a  Versatec  color  plotter  available  for  use  to  this  project.  The  department 
has  a  CONVEX  Cl  supercomputer,  which  is  part  of  the  Stanford  Exploration  Project  (SEP). 
Finally,  we  have  direct  access  to  three  CRAY  computers:  One  through  Stanford  university  to 
the  San  Diego  supercomputing  center;  another  through  our  DOE  grant  to  Lawrence  Berkeley 
Laboratory’s  CRAY;  and  the  third  through  Chevron  to  their  supercomputing  facility  at  LaHabra. 

We  propose  to  use  the  Ardent  for  the  main  supercoiiputing  platform  in  this  study.  All  of 
the  codes  described  above  are  implemented  on  this  machine.  With  floating  point  speed  of 
approximately  one  fifth  the  Cray,  the  Ardent  provides  much  better  tum-around  and  3D  graphics 
for  easy  analysis  an  interpretation  of  results.  The  Ardent  costs  approximately  $50/hr  for  the 
complete  resource  including,  CPU  time,  disk  storage,  operations,  etc.  or  approximately  1/40  the 
cost  of  remote  Cray  time. 

7.2  Personnel  and  Resumes 

Dr.  Jerry  Harris  (see  enclosed  CV)  received  the  Ph.D.  from  the  California  Institute  of 
Technology  in  1980  in  electrical  engineering.  He  is  an  expert  in  the  area  of  signal  processing 
based  on  wave  physics  and  specializes  in  single-well  and  cross-well  seismic  imaging.  During 
the  past  4  years,  he  was  development  manager  of  British  Petroleum's  cross-well  seismology 
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program  and  was  responsible  for  its  advance  over  all  other  such  programs  in  industry.  He 
joined  the  Department  of  Geophysics  as  Associate  Professor  in  the  fall  of  1988. 

Dr.  Amos  Nur  (see  enclosed  CV)  received  the  Ph.D.  in  geophysics  from  the  Mass.  Inst, 
of  Technology  in  1969.  He  is  an  expert  in  rockphysics  and  rock  mechanics  and  specializes  in 
the  relation  between  seismic  and  reservoir  properties  of  rocks.  He  has  published  over  150 
refereed  papers,  resulting  from  funding  over  the  past  15  years  by  the  Department  of  Energy, 
NSF,  GRI  and  for  Ae  last  1 1  years  by  a  consortium  of  oil  and  oil  field  service  companies.  He 
is  professor  of  geophysics  and  also  chairman  of  the  department. 
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7.3  CVs  and  Bibliographies 
Jerry  M.  Harris 

Associate  Professor  of  Geophysics,  Stanford  University 


Year  and  Place  of  Birth: 

1951;  Sardis,  Mississippi 

Education: 

B.S.,  1973;  Electrical  Engineering,  University  of  Mississippi 
M.S.,  1974;  Electrical  Engineering,  California  Institute  of  Technology 
Ph.D.,  1980;  Electrical  Engineering,  California  Institute  of  Technology 

Professional  Positions: 

Staff  Engineer,  Communications  Satellite  Corporation;  1974-77. 

Research  Specialist,  Exxon  Production  Research  Company;  1980-84. 

Staff  Geophysicist,  Standard  Oil  Production  Company;  1984-1988. 

Associate  Professor  of  Geophysics,  Stanford  University;  1988  -  present. 

Professional  Activity: 

Member,  Society  of  Exploration  Geophysics 
Member,  Institute  of  Electrical  &  Electronics  Engineers 
Member,  Society  of  Petroleum  Engineers 

U.S.  Delegate,  International  Consultative  Committee  on  Radio  (CCIR),  Geneva;  1976. 

Research  Interests  and  Activities: 

Seismic  and  electromagnetic  wave  propagation 
Adaptive  high  resolution  signal  processing 
Wave  equation  imaging  and  inversion 
High  resolution  experimental  techniques 

Honors  and  Awards: 

Coming  Fellow,  CalTech;  1974. 

Academic  Reward  for  College  Scientists,  CalTech;  1978. 

Hughes  Fellow,  Caltech,  1979-80. 

Member,  Engineering  honor  society.  Eta  Kappa  Nu  &  Tau  Beta  Pi. 

Outstanding  Graduate  Achievement  Award,  University  of  Mississippi;  1982. 


19. 


J.  M.  Harris'  Recent  Publications 


D.  J.  Fang  and  J.  M.  Harris,  1979.  Precipitation  Attenuation  Studies  Based  on 
Measurements  of  the  ATS-6  20/30  GHz  Beacons  at  Clarksburg,  MD,  IEEE  Trans. 
Antennas  Prop.,  vol.  AP-27,  January  1979,  pp  1-11. 

W.  Wells  and  J.  M.  Harris,  1981.  Multiple  Scattering  of  Collimated  Irradiance,  J. 
Optical  Society  of  America,  vol.  71,  no.  3,  March  1981,  pp  243-249. 

J.  M.  Harris,  1987.  Diffraction  tomography  with  discrete  arrays  of  sources  and 
receivers,IEEE  Trans.  Geoscience  and  Remote  Sensing,  Vol.  GE-25,  No.  4,  pp.  448-455. 

G.  A.  McMechan,  J.  M.  Harris,  and  L.  Anderson,  1987.  Cross-well  tomography  for 
strongly  variable  media  with  application  to  scale  model  data.  Bulletin, 

Seismological  Society  of  America,  Vol.  77,  No.  6,  pp.  1945-1960. 

J.  M.  Harris  and  A.  G.  Ramm,  1988.  Two  Dimensional  Inversion  of  Well-to-Well  Data, 
Appl.  Math.  Lett.,  Vol.  1,  No.  2,  pp.  127-131. 

L.  Z.  Hu,  G.  A.  McMechan,  and  J.  M.  Harris,  1988.  Acoustical  pre-stack  migration  of 
cross-well  data.  Geophysics,  Vol.  53,  No.  8,  pp.  1015-1023. 

J.  M.  Harris,  1989.  Point-to-Point  Bom  inversion  of  spherical  wave  data,  to  appear 
in  Geophysics. 
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VITAE:  AMOS  M.NUR 


Amos  M.  Nur,  Bom  Febmary  9, 1938. 

Position: 

Professor  of  Geophysics,  Chairman  of  the  Department  of  Geophysics. 

Education: 

B.  Sc.,  Geology,  Hebrew  University 

Ph.D.,  Geophysics,  M.I.T.  Cambridge,  Mass.  1969. 

Professional  Positions: 

Research  Associate  in  Geophysics,  MIT.  1969-70. 

Assistant  Professor  of  Geophysics,  Stanford  University,  1970-74. 

Associate  Professor  of  Geophysics,  Stanford  University,  1974-79. 

Visiting  Professor,  Weizmaim  Institute  of  Science,  1974-77. 

Professor  of  Geophysics,  Stanford  University,  1979- 
Chairman,  Department  of  Geophysics,  Stanford  University,  1986- 

Professional  Activity: 

Member,  Lunar  Sample  Review  Board,  1972-73. 

Associate  Editor,  Journal  of  Geophysical  Research,  1974-77. 

Member,  Committee  on  Seismology,  National  Research  Council,  1974-77. 

President,  Tectonophysics  Section,  AGU,  1976-77. 

Member,  National  Science  Foundation  Earth  Sciences  Advisory  Board,  1974-77. 
Member,  National  Academy  of  Science  Earthquake  Delegation  to  the  Peoples  Repu¬ 
blic  of  China,  1976. 

Delegate,  Japan-US.  Earthquake  Prediction  Symposium,  Tokyo,  1977. 

Director,  Stanford  Rockphysics  and  Borehole  Geophysics  Project  (SRB)  1977- 
Associate  Editor,  Tectonophysics,  1979-1987. 

Member  and  Chairman,  Earth  Sciences  Review  Panel,  LBL,  1980-81. 

Member,  Day  Medal  Committee,  G.S.A,  1981. 

Member,  Bucher  Medal  Committee,  A.G.U.  1981. 

Member,  Hydrocarbon  Research  Drilling  Committee,  NRC,  1987-88. 

Invited  speaker  at  1988  SEG  Aimual  Meeting  Research  Symposium. 

Member,  AGU  Macelwane  award  committee,  1988. 

Member,  advisory  board.  Department  of  Earth  Science,  Purdue  University,  1988. 
Honors  and  Awards: 

Research  Fellow,  A.P.  Sloan  Foundation,  1972-74. 

J.  B.  Macelwane  Award,  American  Geophysical  Union,  1974. 

Mellon  Fellowship,  Stanford  University,  1974-75. 

Newcomb  Cleveland  Prize,  A.A.A.S.  1975. 

Fellow,  American  Geophysical  Union,  1976. 

Fellow,  Geological  Society  of  America,  1980. 

University  Fellow,  Stanford  University,  1980. 

Fellow,  Sackler  Institute,  Tel  Aviv  University,  1985- 
Wayne  Loel  Professorship  of  Earth  Sciences,  1988. 

Invited  Snider  Lecturer,  Erindale  College,  University  of  Toronto,  Canada,  1988. 
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Bicentennial  Fellow,  Universidad  de  los  Andes,  Merida,  Venezuela,  1988. 
Visiting  Professor,  University  of  Tokyo  and  ERI,  1988. 

Professional  Societies: 

American  Geophysical  Union. 

Geological  Society  of  America. 

Society  of  E:q)loration  Geophysicists. 

Seismological  Society  of  America. 

American  Association  for  the  Advancement  of  Science. 

Society  of  Petroleum  Engineers. 

New  York  Academy  of  Sciences. 


PUBUCATIONS:  AMOS  NUR 


1968.  Gene  Simmons  and  Amos  Nur.  GRANITES:  RELATION  OF  PROPERTIES  IN 
SITU  TO  LABORATORY  MEASUREMENTS.  15  November.  V162,  pp  789-791. 

1969.  Amos  Nur.  THE  EFFECT  OF  FLUID-FILLED  PORES  ON  SEISMIC  VELOCI¬ 
TIES  OF  ROCKS.  Earthquake  Notes.  VXl,  pp  1-24. 

1969.  Gene  Simmons  and  Amos  Nur.  GRANITIC  ROCK:  PROPERTIES  IN  SITU. 
Science,  V165.  pp  203. 

1969.  Amos  Nur  and  Gene  Simmons.  THE  EFFECT  OF  SATURATION  ON  VELOaTY 
IN  LOW  POROSITY  ROCKS.  Earth  and  Planetary  Science  Letters  7,  pp  183-193. 

1969.  Amos  Nur  and  Gene  Simmons.  THE  EFFECT  OF  VISCOSITY  OF  A  FLUID 
PHASE  ON  VELOCITY  IN  LOW  POROSITY  ROCKS.  Earth  and  Planetary  Science 
Letters  7,  pp  99-108. 

1969.  Amos  Nur  and  Gene  Simmons.  STRESS-INDUCED  VELOCITY  ANISOTROPY  IN 
ROCK:  AN  EXPERIMENTAL  STUDY.  JGR.  V74,  #27,  pp  6667-6674.  December  15. 

1970.  A.  Nur  and  G.  Simmons.  THE  ORIGIN  OF  SMALL  CRACKS  IN  IGNEOUS 
ROCKS.  Int.  J.  Rock  Mech.  Min.Sci.  V7,  pp  307-314. 

1970.  Hiroo  Kanamori,  Amos  Nur,  D.  Chung,  D.  Wones  and  Gene  Simmons.  ELASTIC 
WAVE  VELOCITIES  OF  LUNAR  SAMPLES  AT  HIGH  PRESSURES  AND 
THEIR  GEOPHYSICAL  IMPLICATIONS.  Science.  30  January,  V167.  pp  726-728. 

1970.  Ki-iti  Harai  and  Amos  Nur.  RELATIONSHIP  AMONG  TERRESTRIAL  HEAT 
FLOW,  THERMAL  CONDUCTIVITY,  AND  GEOTHERMAL  GRADIENT.  JGR. 
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